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ABSTRACT 

We analyze the observed distribution of the orbital eccentricity and period of 
binary radio pulsars in globular clusters using computational tools to simulate 
binary-single star interactions. Globular clusters have different groups of pulsars 
arising from separate interaction scenarios. Intermediate eccentricities of cluster 
pulsars can be mostly accounted by fly-bys although locally lower stellar densities 
at pulsar positions may alter the situation. Very high eccentricities are likely to 
be results of exchanges and/or mergers of single stars with the binary companion 
of the pulsar. 

Subject headings: pulsars: general — globular clusters: general 

1. Introduction 

Sustained high sensitivity search of globular clusters (GC) for radio pulsars with im- 
proved pulsar-search algorithms have yielded so far 140 pulsars in 26 GCd3 (Camilo & Rasio 
2005 and references therein). Radio pulsars can be timed more easily and accurately by 
ground-based telescopes over well-separated epochs, as the underlying neutron stars (NS) 
are less prone to noise from episodically varying accretion torques (as for X-ray pulsars). 
These lead to the easier measurement of orbital parameters. Binary systems can provide 
an important source of energy for GCs, since the binding e nergies of a few , very close bi- 



naries can approach that of a moderately massive host GC (IHut et al. 20031 ) and therefore 
can have dynamical effects on the cluster's evolution. On the other hand, the observable 
parameters of the neutron stars and their binary companions in the clusters, such as spin, 
orbital period and eccentricity, projected radial position in the cluster, companion mass and 



information on these pulsars are found from P. Freire's webpage updated August 2008, 



http://www.naic.edu/~pfreire/GCpsr.html, compiled from radio timing observations by many groups. 
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their distributions provide a valuable test-bed to examine the theoretical scenarios of forma- 
tion and evolution of recycled pulsars. These parameters can provide a tracer of the past 
history of dynamical interactions of the binary NSs in individual GCs. We discuss below 
the distribution of orbital eccentricities and periods of GC pulsars in the light of interaction 
of stars with binaries already formed or in the process of formation inside globular clusters, 
due to the high stellar density in their cores. The observed distributions are examined with 
analytical results and our numerical experiments of scattering of stars simulated by direct 
N-body integration tools. 

2. Observed data 

Although there are a number of low mass x-ray binaries (LMXBs) in the GCs, their 
orbital period and eccentricity information is poorly known compared to radio pulsars; more- 
over the initial eccentricity of formation is quickly quenched due to tidal coupling with the 
Roche lobe filling companion. Hence we concentrate entirely on binary radio pulsars whose 
companions are usually smaller than their Roche lobes, and less well coupled tidally. Of 
the 140 radio pulsars in the GCs, 74 binaries (a large fraction of which are millisecond 
pulsars), 59 are isolated and 7 others have no published timing/orbital solutions. Orbital 
parameters are not well determined for one binary (PSR J2140— 2310B in M30), only the 
lower limit of P or b and e are known as P or & > 0.8 days and e > 0.52 and thus excluded 
from the present study. In Fig [TJ we plot e Vs P or b for 73 binaries in GCs with known 
orbital solutions. Logarithmic scales in both eccentricity and orbital periods are chosen, as 
the enormous range of both variables and the regions occupied by observed pulsars are less 
obvious in linear scales. Observed binary radio pulsars in GCs can be categorized into three 
groups: I) 21 pulsars with large eccentricity (1 > e > 0.01); II) 20 pulsars with moderate 
eccentricity (0.01 > e > 2 x 10~ 6 ); III) 32 pulsars with small eccentricity (e ~ 0). Using the 
'R' statistical package Efk-means test, we show that these groups are distinct at a statistically 
significant leve|E 

In the database, several pulsars' orbital eccentricities have been listed as zero, but we 
assign them an arbitrarily small value of e = 3 x 10~ 7 . Note that the smallest eccentricity 



2 www . R-pro j ect .org 

3 Here we assign random eccentricities of the group III pulsars in such a way that they remain below the 
"limit of timing sensitivity" line. We got 3 clusters of sizes 18 (medium eccentricity), 34 (low eccentricity), 
21 (high eccentricity) with the sum of squares from points to the assigned cluster centers as 10.2, 13.2 and 
19.9 respectively whereas the squares of inter-cluster distances are : d\ 2 = 7.6, c?23 = 30.5, d\ x = 10.8. 
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measurable is determined by how well the or bit of the binar y is sampled and the overall timing 
accuracy achieved for radio pulse arrivals (jPhinney 1992T ) . The timing accuracy translates 
into an upper limit on the smallness of the eccentricity This is displayed as the brown solid 
line in the lower left corner of Fig [T] following the functional relation: e m i n = (St) / (asin i/c) = 



4TT 2 cSt/ 



sin i 



m p + m c 



,1/3 



2/3 
orb 



Here i is the inclination angle of the binary, m p is the 
pulsar mass, m c is the companion mass and St is the timing accuracy. We take m p = 1.4 M , 
m c = 0.35 M Q , i = 60°, St — 1 /i.sec for the "limit of timing sensitivity" line. 

Obse rvational selection effe cts may be influencing the distribution of GC pulsars seen 
in Fig [T] (jCamilo fc Rasio 20051 ). the most important selection effect operates towards the 
left of the diagram: it is more difficult to detect pulsars with larger DM and/or shorter spin 
periods, especially millisecond pulsars in short orbital period and highly eccentric binaries. 
Another important selection effect is due to distance, since only the brightest pulsars can be 
observed at large distance. 



3. Two- and three-body stellar interactions in globular clusters 



The presence of a large number of L MXBs in GCs compared to the galactic field had 
led to the suggestion (IFabian et al. 19751 ) that a binary is formed by tidal capture of a non- 
compact star by a neutron star in the dense stellar environment of the GC cores. If stable 
transfer of mass and angular momentum ensued from the comp anion star, this co uld lead to 
recycled pulsars in binaries or as single millisecond pulsars e.g. (jAlpar et al. 19821 ). However 
tidal capture of a neutron star by a low mass main sequence star can lead to large energies 
being deposited in tides. The resultant structural readjustments of the star in response 
to the dissipation of the modes could be very significant in stars with either convective or 
radiative damping zones (IRay et al. 19871 : iMcMillan et al. 19871 ) and the companion star 
can undergo a size "inflation" due to its high tidal luminosity which may be much larger 
than that induced by nuclear reactions in the core. Efficiency of viscous dissipation and 
orbit evolution is crucial to the subsequent evolution of the system as viscosity regulates the 
growth of oscillations and also the extent the extended star is bloated and shed. A significant 
fraction of the encounters lead to binaries that either become unbound as a result of de- 
excitation or heating from other stars in the vicinity or they are scattered into orbits with 
large pericenters (compar ed to the size of the non-compact star) due to angular momentum 
transfer from other stars (IKochanek 19921 ). For a recent summary of the formation channels 
of retained neutron star binaries in GCs obtained via population synthesis, see Ivanova et al. 
(2008). 



If the local binary fraction is substantial, the single star - binary interaction can ex- 
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ceed the encounter rate between single stars by a large factor (jSigurdsson fc Phinnev 19 93). 



The existence of a significant population of primordial binaries in GCs ( jYan &: Mateo 1994 



Pryor et al. 19891 ) indicate that three body processes have to be accounted for in any dynam- 
ical study of binaries involving compact stars. For a literature summary of the constraints 
on the binary fraction in globular clusters see Davis et al. 2008. 

An encounter between a field star and a binary may lead to a change of state of the 
latter, e.g.: i) the original binary may undergo a change of eccentricity and orbital period 
but otherwise remain intact - a "fly-by" interaction; ii) a member of the binary may be 
exchanged with the incoming field star, forming a new binary - an "exchange" process; iii) 
two of the stars may collide and merge into a single object, and may or may not remain 
bound to the third star - a "merger" process; or iv) all three stars become unbound - an 
"ionization" process. 



4. Fly-by, exchange- and merger- collision induced eccentricity 

Formation scenario of millisecond pulsars in primordial binaries suggests that their 
eccentricity should be very small, ~ 10~ 6 — 10~ 3 (Phinney 1992). But inside globular clusters, 
they may acquire eccentricity through interactions with single stars. Rasio & Heggie (1995), 
Heggie & Rasio (1996) studied the change of orbital eccentricity (Se) of an initially circular 
binary following a distant encounter with a third star in a parabolic orbit. They used 
the secular perturbation theory, i.e., averaging over the orbital motion of the binary for 
sufficiently large values of the peri centre distance r p , where the encounter is quasi-adiabatic 
and used non secular perturbation theory for smaller values of r p where the encounter is 
non-adiabatic. In the first case be varies as a power law with r p /a and in the second case 
be varies exponentially with r p /a (a is the semi- major axis of the binary). The power law 
dominates for e < 0.01 and the exponential dominates for e > 0.01. They estimated the 
cross-sections (a) for eccentricity changes and calculated time-scales for eccentricity changes 
as t = 1/ {rate) = 1/ < nav > where n is the number density of the stars and v is the 
velocity of the incoming star. The expressions of the timescales for fly-by are (see Rasio & 
Heggie 1995): 

= 4 x 10 n n^ l v 10 p- 2 b /3 e 2/5 for e < 0.01 

t fly = 2 x 10 u ri4 ^ioP'T [- ln(e/4)]" 2/3 for e > 0.01 

where is the number density (n) of single stars in units of 10 4 pc -3 and viq is the velocity 
dispersion (i>) in units of 10 km/sec in GCs; P or b is the orbital period in days giving tfi y in 
years. 
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The eccentricities of the binary pulsars in GCs are likely to be due to binary single 
star interactions when the interaction timescale is less than the binary age. We take the 
maximum age of the binaries in a GC to be the globular cluster ages which are ~ 10 10 
years. The value of thoAu whic h determines t/i y varies from 0.0024 - 2.167 for different 



globular clusters (IWebbink 19851 ). We grouped them according to the values of f 10/^4 and 
calculated tfi y with the mean values of v 10/^4 for each of six groups. In Fig. [TJ we plot the 
isochrones of fly-by encounters (tfi y ) in the e — P or b plane for all six groups using the above 
expressions. Pulsars which lie inside or outside the globular cluster core in the projected 
image are depicted with different symbols and the colors are the same as that used to draw 
corresponding tfi y contours. If tfi y > 10 10 years for a particular binary, then it would not 
have interacted and it would preserve its original eccentricity. If tfi y < 10 10 years for a 
particular binary, then it could be eccentric due to fly-by interactions which is the case for 
most of the eccentric GC binaries (Fig. [TJ. So, many GC pulsars' orbital eccentricities can be 
explained by fly-bys. However the local stellar densities at the pulsar positions (especially 
with positional offset from cluster cores) may indicate a higher vio/n^ value from central 
values used to calculate the isochrones. Thus these pulsars might be in regions where the 
effective timescale for fly-by induced eccentricity is larger than the Hubble time. Moreover, 
for pulsars with high eccentricities (e > 0.1, a majority of the group I pulsars) a very close fly- 
by or multiple fly-bys are necessary if the initial binary was circular (see Heggie & Rasio 1996 
and Camilo & Rasio 2005 for a discussion). In these cases, other processes, e.g. exchange and 
merger may produce these systems more naturally. We discuss this below with numerical 
simulations. 

We have used the STARLAE0 software (which reports the results for both resonant and 
non-resonant encounters for exchange and merger processes) to perform numerical simula- 
tions. We discuss the results for pulsar binaries found in Ter 5, since it has the lowest value 
of the parameter v 10/^4 where all classes of encounters can take place. We find that the 
"phase-space" of the encounter reactions of different kind in 47 Tuc and other clusters, is 
not as large as that of Ter 5. In Fig. [2J we plot P or b,in of the initial binary along the top 
x-axis and P or bjin along the bottom x-axis. Porbjin is obtained from P or b,in putting A = in 
the relation aji n = [a,i n m a mb/ 'm^ (1 — A)] where mi and m-i are masses of the members 
of the initial binary, m 3 is the mass of the incoming star, m a and m& are masses of the mem- 
bers of the final binary, A is the fractional change of binary binding energy. For exchange, 
m a = mi, mb = for merger, m a = mi, nib — m 2 + WI3. 

It is clear from the scatter plots (Fig. [2D that the final binaries will most probably have 
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e > 0.1 if they undergo either exchange or merger events. Six high eccentricity (e > 0.1) 
binaries in Terzan 5 are also shown (red in colors, symbols same as used in Fig[T]) in this plot. 
All of them might result from exchange interactions with either a normal mass companion 
(~ 0.40 M Q ) or with a low mass companion (~ 0.16 M ). PSR U, X and Z might even come 
from exchange with ultra low mass companion (~ 0.024 M ). Note that all the exchange 
conclusions here are really guidelines, since they are based on stars of a single mass (0.33 M & ) 
exchanging into the systems, and some of the minimum inferred companion masses are quite 
different from 0.33 M Q . Mergers with 0.40 M & (initial) companions and incoming 0.33 M & 
stars are problematic because of the high final companion masses. Q is the only system for 
which this might be possible, and that would require a fairly small orbital inclination angle. 
Similar problems apply to initial companion masses of 0.16 M Q for all except Q and U. 
Finally, while the mass restriction is not really a problem for any of the lower-mass systems 
in the ultra-low-mass case, the small orbital periods imply that too long a time would have 
to pass before a suitable encounter took place. 

In Fig [U there is a cluster of three pulsars with 0.01 < e < 0.1 and 60 < P or b < 
256 d - NGC1851A, M3(D) and Ter 5(E), all with companions of mass in the range m c = 
0.21 — 0.35 M Q assuming % = 60° (except NGC 1851A which has a 1.12 M companion). 
These are possibly w hite dwarf cores of red giant companions that overflowed Roche lobe 



( IWebbink et al. 19871 ). Such binaries would normally have the "relic" eccentricities ~ 10~ 4 
(Phinney 1992). The above binary pulsars with their presently mildly high eccentricities, 
have undergone fly-by encounters with field stars, rather than exchange reactions, which 
would produce very high eccentricities e > 0.1. Ter 5 E lies outside the core but even 
then, the density may have been high enough to allow a strong fly-by interaction. It could 
also have been ejected out of the high density core after a strong interaction. In addition, 
the pulsar B1620— 26 (in M4) occupies approximately the same region of the phase space. 
However, this system is most likely a triplet system with a planet sized third body and their 
interactions lead to the characteristics of the inner orbit (Thorsett et al. 1999, Ford et al. 
2000, Sigurdsson et al. 2003 and references there in). 

Another set of three millisecond pulsars have 0.01 < e < 0.1, 2 < P orb < 10 d; Ter 5 
(W), 47 Tuc (H) and NGC6440 (F). These clusters have low values of v xo/n^, and so fly-by 
encounters in these clusters would be efficient and could generate these eccentricities in GCs, 
even if their progenitor binaries had short orbital periods and had sub-giant companions of 
the NSs. Alternately, these binaries could also have been formed by fly-by interactions from 
a presently less abundant longer period 2 < P orb < 10 d cluster of "intermediate eccentricity" 
binaries to the right of the pulsars seen in the middle of Fig [TJ 

The "intermediate eccentricity" binaries (group II : 0.01 > e > 2x 10 -6 ), could have been 
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generated by fly- by encounters with low (or "zero" ) eccentricity progenitor pulsars below the 
line of "timing sensitivity limit" (group III pulsars). Some of the shorter P or b binaries would 
again be circularised by gravitational radiation (see t gr contours in Fig [T] calculated using 
the formalism of Peters & Mathews 1963). The progenitor group III pulsars, themselves 
occur in regions of favorable fly-by encounters inducing higher eccentricities. These nearly 
circular binaries have their P or b in the range of 0.06 — 4 days among which the tighter 
binaries can not be progenitors of exchanges/mergers whenever interaction timescales are 
greater than 10 10 years (see Fig. [2]). The minimum P^b, i n above which they can undergo 
exchange/merger interactions decreases as m-i decreases. About the origin of the group III 
pulsars themselves, we note that Camilo & Rasio (2005) discuss the dynamical formation of 
ultra-compact binaries involving intermediate mass main sequence stars in the early life of 
the GC. These companions must have been massive enough (beyond the present day cluster 
turn-off mass of 0.8 M Q ) so that the initial mass transfer became dynamically unstable, 
leading to tight NS-WD binaries through common envelope evolution. Alternately, present 
day redgiant and NS collisions lead to a prompt disruption of the redgiant envelope and the 
system ends up as eccentric NS-WD binary (Rasio & Shapiro 1991). NS-WD binaries can 
be circularised to group III by gravitational wave radiation if P orb < 0.2 days (see Fig. [TJ 
here and also discussions by Camilo & Rasio 2005). 

In conclusion, we find that the presently observed orbital eccentricity and period data 
of GC binary pulsars are largely consistent with numerical scattering experiments on stellar 
interaction scenarios of fly-bys, exchanges and mergers with typical field stars characterizes 
by the central regions of the GCs. The efficiency of fly-by interaction is subject to the 
local stellar density at the location of the initially circular binaries. Exchange and merger 
interactions induce the highest range of eccentricities (1 > e > 0.1) and may be operative in 
different GCs. 

We thank Roger Blandford, Avinash Deshpande and Shri Kulkarni for discussions, Sayan 
Chakraborti for comments on the manuscript and the anonymous referee for constructive 
suggestions. We thank the STARLAB development group for the software and the ATNF 
pulsar group and Paulo Freire for pulsar data bases. This research is a part of 11th plan 
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Fig. 1. — Cluster binary pulsars in the e — P or b plane with contours of tfi y = 10 10 yrs (solid 
lines) and tfi y = 10 s yrs (dashed lines) for different values of fio/^4- Contours of t gr = 10 10 
yrs (solid lines) and t gr = 10 8 yrs (dashed lines) for binaries with m p = 1.4 M & and 
m c = 0.35 M or 0.16 M was calculated using the formalism of Peters & Mathews (1963). 
Pulsars with projected positions inside the cluster core are marked with +, those outside 
the cluster core with x and the pulsars with unknown positions with o (see footnote 1 and 
S. Ransom's webpage at www.cv.nrao.edu/~sransom/). The color scheme for tfi y and t gr 
contours are shown on the right. Solid red curve for tfi y = 10 10 years for Vy^jn^ = 0.0024 
(Terzan 5) is outside the range plotted. Individual pulsars are marked with same colors as 
v 10/714 values of their host GCs. A pulsar is located on the upper left half of the corresponding 
tfiy = 10 10 years line (e.g. the pulsar shown by a black diamond in the middle of the figure) is 
most likely a primordial binary unless it is in the range of eccentricities typical for exchanges 
and merger. Six high eccentricity (e > 0.1) binaries in Terzan 5 are marked with their names. 
Different groups of GCs according to the increasing values of v 10/^4 are : (i) Terzan 5 (red) 
(ii) NGC 6440, M30, NGC 1851, M62, M15 (green) (iii) NGC 6441, NGC 6544, 47Tuc (blue) 
(iv) M28, NGC 6342, NGC 6752, NGC 6760, NGC 6539, NGC6397 (magenta) (v) M4, M5, 
M3, M22 (cyan) (vi) M71, M13, NGC 6749, M53 (black). 
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Fig. 2. — Time scales (denoted by '+') and final eccentricity distributions (scatter-plot of 
points) with initial and final orbital periods (A = 0) for exchange (purple points on the 
left panel) and merger (green points on the right panel) interactions with different stellar 
parameters. Each set of panels had 5000 trial densities for each orbital period for the mass 
combination shown (typiclly 15570 total scatterings led to 10960 fly-bys, 3722 exchanges, 
982 two mergers and 6 three- mergers for each P or b,in in set 2). We plot P or b,in along the top 
x-axis and P or bjin along the bottom x-axis. The left y axis gives the final eccentricities while 
the right y axis gives the time scales of interactions. Vertical orange lines give boundaries 
of orbital periods where interaction time scales < 10 10 yrs. Six high eccentricity (e > 0.1) 
binaries in Terzan 5 are also shown (colored red, symbols same as used in Fig [1] regarding 
their positions) in each case. 



